Enhancement of the retrapping current of superconducting microbridges of finite 

length 
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We theoretically find that the resistance of a superconducting microbridge/nanowire decreases 
while the retrapping current I r for the transition to the superconducting state increases when one 
suppresses the magnitude of the order parameter |A| in the attached superconducting leads. This 
effect is a consequence of the increased energy interval for diffusion of the 'hot' nonequilibrium 
quasiparticles (induced by the oscillations of |A| in the center of the microbridge) to the leads. 
The effect is absent in short microbridges (with length less than the coherence length) and it is 
relatively weak in long microbridges (with length larger than the inelastic relaxation length of the 
nonequilibrium distribution function). A nonmonotonous dependence of I r on the length of the 
microbridge is predicted. Our results are important for the explanation of the enhancement of 
the critical current and the appearance of negative magnetoresistance observed in many recent 
experiments on superconducting microbridges/nanowires. 

PACS numbers: 74.25.Op, 74.20.De, 73.23.-b 



I. INTRODUCTION 

Recently, several experimental groups observed a 
negative magnetoresistance (NMR) of superconducting 
nanowires/microbridges at temperatures lower than the 
critical temperature [lKD]. In Refs. [3|-l6j authors demon- 
strated that in their case the effect is connected with the 
suppression of superconductivity in the bulk supercon- 
ductor caused by the applied magnetic field. Moreover 
in Refs. 0-Q it was shown that this NMR appears to- 
gether with an enhancement of the critical current of the 
nanowires in weak magnetic fields. 

Presently there exist several theories [||-[l2[ which pro- 
pose different mechanisms for the NMR. Ref. Q claims 
that it is connected with the suppression of a new chan- 
nel of dissipation in superconducting wires by an ap- 
plied magnetic field, while Refs. [3, [l(| argue that 
the suppression of the intrinsic pair-breaking resulting 
from total 'spin-flip'+'non-spin-flip' rate is responsible 
for the effect. In Ref. [ll| the stabilization of the su- 
perconducting phase due to magnetic field induced nor- 
mal metal/superconductor boundaries at the ends of the 
microbridge was offered as the main mechanism for the 
negative magnetoresistance while in Refs. [T], UH, [l3[ the 
effect was explained as due to a decrease of the charge 
imbalance decay length Xq in weak magnetic fields. 

Taking into account the strong dependence of both the 
NMR and the enhancement of I c on the length of the su- 
perconducting microbridge/nanowire (i.e. the effect does 
not exist in relatively long and short samples - see Refs. 
0-Q) we argue that the proposed mechanisms as put 
forward in Refs. [8l-fll| are not relevant for these exper- 
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iments. Indeed their applicability is not limited by the 
length of the superconductors. Due to the same reason 
a decrease of Xq [l|, UH, EH cannot explain the length 
dependence of the effect because it should lead to NMR 
in samples of arbitrary length. Besides at currents close 
to the depairing current , Xq starts to depend on the cur- 
rent [TH and it smears its dependence on H. As a result 
this mechanism cannot explain the increase of I c in weak 
magnetic fields. 

In Ref. [IH another mechanism for the enhancement of 
I c was proposed connected with the complete suppression 
of the order parameter in the superconducting leads. But 
this mechanism is not able to explain the monotonous 
enhancement of I c in weak magnetic fields (see [2, la, lf| ) 
when the superconducting leads are still in the supercon- 
ducting state with weakly suppressed order parameter. 

Here, we offer a new mechanism that leads to the neg- 
ative magnetoresistance of microbridges/nanowires and 
the enhancement of the critical (retrapping) current. Our 
explanation is based on the following idea: even a weak 
suppression of the order parameter in the leads opens 
new energy channels for the diffusion of the 'hot' quasi- 
particles from the microbridge where they are induced 
by oscillations of the order parameter [15|]. As a result 
the effective 'temperature' of the quasiparticles in the 
microbridge decreases and the retrapping current for the 
transition to the superconducting state increases. Simul- 
taneously it leads to decrease of the resistance of the 
microbridge at fixed current. Our proposed mechanism 
leads to no effect in very short and very long micro- 
bridges/nanowires (in agreement with the experiments 
Besides we find that the retrapping current is a 
nonmonotonous function of the length of the microbridge. 

The paper is organized as follows. In section II we 
present the theoretical model. The current- voltage (IV) 
characteristics of microbridges of different length are pre- 
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sented in Section III for different values of the order pa- Ojj 
rameter in the leads. In section IV we discuss our result 
and in Sec. V we present our conclusions. 
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FIG. 1: (Color online) (a) Schematic illustration of the fi- 
nite width Dayem microbridge connected with wide super- 
conducting leads, (b) Dependence of the magnitude of the 
order parameter along the center line of the microbridge for 
different magnetic fields, (c) Contour plot of |A| in the leads 
and microbridge at different H (blue(red) color corresponds 
to low(high value for |A|)). The length of the full system 
(leads+microbridge) is 27£, width 15£, length of microbridge 
is 7£ and width of microbridge is £. 



To study the resistive state of the microbridge we use 
the set of equations (time-dependent Ginzburg-Landau 
equation for the superconducting order parameter cou- 
pled with the kinetic equations for the quasiparticle 
distribution function and the Usadel equations for the 
Green's functions) derived for 'dirty' superconductors 
near the critical temperature of the superconductor [16l — 
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Here A = |A|e^ is a complex order parameter, £g L = 
TrhD/8k B T c and A 2 GL = 8ir 2 (k B T c ) 2 /7C(3) are the zero 
temperature Ginzburg-Landau coherence length and the 
order parameter correspondingly. Q — (d<fi/dx — 2eA/hc) 
is a quantity which is proportional to the superfluid ve- 
locity (v s = DQ) with A the vector potential taken 
in the Landau gauge, ip is the electrostatic potential, 
/i( e ) = fl( e ) + ^/L( e ) is tne longitudinal and / T (e) 
is the transverse parts of the quasiparticle distribution 
function 2/(e) = 1 — /l(c) — /r(e) (in equilibrium Jl = 
fl(e) = tanh(e/2fc B T 1 ), /£(e) = 0). Ni, N 2 , R 2 are the 
spectral functions which are determined by the Usadel 
equation for the normal a(e) = cosG = iVi(e) + i-Ri(e) 
and anomalous ft = /3e^, ft = ft~ # (/3(e) = sin 6 = 
N 2 (e) + iR 2 (e)) Green functions. 

Nonequilibrium corrections to the quasiparticle dis- 
tribution function enters Eq. (3) via the potentials 
$i = R 2 Sf L de/\A\ and $ 2 = N 2 f T de/\A\. Eqs. 
1 (a,b) are coupled due to the finite spectral supercurrent 
i £ = i?e(ftVft - ftVft)/2 = 2N 2 R 2 Q. 

The advantage of Eqs. (1-3) in comparison with the 
ordinary or the extended |18j time-dependent Ginzburg- 
Landau equations is that they allow to take into account 
nonlocal effects. Here, under nonlocality we mean that if 
in one place of the superconductor and in some moment 
in time the quasiparticle distribution function /(e) be- 
comes nonequilibrium (due to some kind of perturbation) 
then nonequilibrium correction to /(e) will be nonzero 
over a distance ~ ft„ = (Dr^) 1 / 2 around that point and 
during a time ~ Ti n after turning off the perturbation. 
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Before solving Eqs. (1-3) numerically we calculate the 
order parameter in a model 2D system (see Fig. 1) in 
the stationary state in the presence of an applied mag- 
netic field. This result demonstrates the suppression of 
|A| in the leads by increasing H and the weak influence 
on |A| in the microbridge. It also shows that instead of 
the 2D model we may use a ID model where the sup- 
pression of the order parameter in the superconducting 
leads could be simulated by introducing locally a lower 
critical temperature. Thus we may use the model of Ref. 
(l5j where we introduce a different critical temperature 
in the leads (see Fig. 2). Correlation between T' c and H 
is clear: smaller T' c corresponds to larger H. 
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order parameter is slightly suppressed in the leads. We 
should stress here that the critical current I c of the mi- 
crobridge (at which the superconducting state becomes 
unstable) monotonically decreases with decreasing T' c . 
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FIG. 2: Schematic illustration of the model system. Critical 
temperature in the leads is smaller than the critical temper- 
ature of the microbridge. 

The self-consistent set of Eqs. (1-3) was solved numer- 
ically using the method and boundary conditions pre- 
sented in Ref. [l5|. Further, we use the following di- 
mensionless units. The order parameter is scaled by 
Ao (Ao = 1.76fesTc), distance is in units of the zero 
temperature coherence length £o = ^/HDJAq, time in 
units of io = ^/Ao and temperature in units of the crit- 
ical temperature T c . The current is scaled in units of 
jo = AoOYi/(£oe) an d the electrostatic potential is in 
units of ipo = Ao/e. It is useful to introduce the di- 
mensionless inelastic relaxation time f in = T in /t which 
is the main control parameter in the model described by 
Eqs. (1-3). For example in MgB 2 f in ~ 4 20], in Nb and 
Pb fi n ~ 40, in Sn fi n ~ 70, in Al fi n ~ 3500 and in Zn 

f in ~2-10 4 0. 

III. RESULTS 

In Fig. 3 we present the IV characteristics of the super- 
conducting microbridge with length L = 21£o calculated 
for different values of the order parameter in the leads 
24] (in the inset we show distribution of the time av- 
eraged order parameter in the microbridge and leads at 
/ = 0.75/ c ). Notice that the retrapping current increases 
and the resistance of the microbridge decreases when the 



FIG. 3: (color online). Current voltage characteristics of the 
superconducting microbridge with length L — 21£o — 4.9£(T) 
calculated at T = 0.92T C and f in = 250. The inset shows 
the distribution of the time-averaged order parameter in the 
microbridge at different T' c . Current is normalized in units of 
the critical current of the microbridge with T' c = T c . 

Our result can be explained as due to the enhanced dif- 
fusion of the 'hot' quasiparticles induced by oscillations of 
the order parameter in the center of the microbridge [TsT ] 
when Ai ea d becomes smaller. Indeed, in this case the 
energy barrier connected with the finite energy gap at 
e < |A|/ ea( 2 decreases and nonequilibrium quasiparticles 
in the wider energy interval can diffuse to the leads. Here 
we should remind the reader that the local energy gap in 
the microbridge is smaller than the local value of the 
order parameter due to the finite supercurrent and the 
spatial variation of |A| [3,[2l[. Therefore, 'hot' quasi- 
particles may diffuse in the energy interval e > |A|j e ad 
even when the local |A| in the microbridge is larger than 
|A| Zead (see inset in Fig. 3). 

In order to illustrate the above effect we show in Fig. 4 
the energy dependence of the nonequilibrium correction 
to fi, in the center of the microbridge averaged over one 
period of Josephson oscillations. In the inset to Fig. 4 
we present the spatial dependence of the time-averaged 
potential <&i which corresponds to the effective 'temper- 
ature' Tejf = T + T c $! of quasiparticles in the micro- 
bridge [15]. One can easily see that with decreasing A; ea( j 
the energy interval where the quasiparticles are 'heated' 
(corresponds to a negative sign of 5Jl) decreases and it 
results in drastic changes of <E>i. 

The found effect depends strongly on the length of the 
microbridge. When L < 2£(T) the value of A/ ea ^ has 
the strongest effect on the value of I r (which itself is 
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FIG. 4: (color online). Energy dependence of the time av- 
eraged S/l in the center of the microbridge with length 
L = 21£ ^ 4.9£(T), f in = 250 and T = 0.92T C calculated 
for I = 0.75J C . 



close to the critical current of the microbridge I c - see 
Fig. 5) while the effect of 'heating' is relatively weak 
(see Ref. [22]). With decreasing A/ ea( j both I c and I r 
monotonically decreases and the resistance increases - see 
Fig. 5. The critical length L c for which I r starts to 
increase depends on the inelastic relaxation time - the 
larger fi n the shorter L c . For example ll£o < L c < 15£o 
for fi n — 250 while 17£o < L c < 21£o for fi n = 60 at 
T = 0.92T C . 



the 'hot' quasiparticles occurs mainly in the microbridge 
and hence the effect of diffusion to the leads is rather 
weak. This is illustrated in Fig. 6 where we present 
IV curves for a relatively long microbridge L — 31£o (— 
0.6Li n at fi n = 250). The change in the IV curves is 
much weaker in comparison with the shorter microbridge 
(compare with Fig. 3) and we found that for L = 51£o ~ 
Li n the effect practically disappears for f, n = 250. 




FIG. 6: (color online). Current voltage characteristics of the 
superconducting microbridge with length L — 31£o — 7.3f (T) 
calculated at T = 0.92T C and fj n = 250. The inset shows 
the distribution of the time-averaged order parameter in the 
microbridge at different T' c . Current is normalized in units of 
the critical current of the microbridge with T' c = T c . 
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FIG. 5: (color online). Current voltage characteristics of the 
superconducting microbridge with length L — 7£o — 1.7£(T) 
calculated at T — 0.92T C and n„ = 250. The inset shows 
the distribution of the time-averaged order parameter in the 
microbridge at different T' c . Current is normalized in units of 
the critical current of the microbridge with T c =T C . 

For microbridges with length L 3> the relaxation of 



There is one interesting effect when L < L c . We find a 
decrease of I r and an increase of the resistance at I ~ I r 
but starting from some current I > I r the resistance of 
the microbridge decreases when A; ea( j is suppressed - see 
Fig. 7. At this length there is a competition between 
the influence of the order parameter in the leads and the 
'heating' of the quasiparticles on the IV curves. At low 
currents J ~ I r the value of I r is determined mainly 
by Ai ea d and I r decreases with decreasing Ai ea d s - At 
larger currents the 'heating' of the quasiparticles becomes 
stronger because it increases with decreasing Josephson 
period. As a result decreasing A; ea d weakens the 'heat- 
ing' effects and the voltage at fixed current decreases. 



IV. DISCUSSION 

A typical feature of many experiments about the 
enhancement of the critical current and the negative 
magnetoresistance is presence of finite resistance of the 
nanowires/microbridges even at low temperatures 0-0| 
which implies a strong influence of fluctuations. It can ex- 
plain the absence of the hysteresis of IV curves observed 
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FIG. 7: (color online). Current voltage characteristics of the 
superconducting microbridge with length L = ll£o — 2.6£(T) 
calculated at T = 0.92T C and fi n = 250. The inset shows 
the distribution of the time-averaged order parameter in the 
microbridge at different T' c . Current is normalized in units of 
the critical current of the microbridge with T' c = T c . 



fore, during a finite time min{ri„, Tdiff ~ L 2 /D} after 
the phase slip event the effective 'temperature' of quasi- 
particles will be larger than the bath temperature and 
the probability of the next phase slip becomes higher. It 
creates the condition for phase slip avalanches. By de- 
creasing the order parameter in the leads one increases 
the flux of 'hot' quasiparticles from the microbridge and 
therefore decreases the effective 'temperature' and the 
probability of phase slip avalanches. The effect should 
be strongest in microbridges with L < Li n where such a 
diffusion is the most effective one. In short microbridges 
with L < 2£(T) the suppression of the order parameter 
in the leads suppresses also A in the microbridge (due to 
the proximity effect - see Figs. l(b,c)) and it gives the 
leading contribution to the probability for phase slips - 
i.e. it considerably increases. Therefore one should ob- 
serve positive magnetoresistance in short microbridges. 
In MoGe Lj„ ~ &0nm (when £o — 137im - see Ref. [23|). 
in Nb L in - 170nm (£o=28 nm 2]) in Pb L m ~ 240nm 
(£o=40 nm), in Al Li n ~ 8/jto (^q—IAQ nm) and in Zn 
Li n — 25/zm (£o=250 nm). It correlates with the length 
of the nanowires/microbridges which were studied in 1- 
7] and where the negative magnetoresistance at low cur- 
rents were observed both near and far from T c . 



in Refs. 0-0] ■ Indeed, it is known for example from 
the theory of Josephson junctions that fluctuations may 
completely destroy the hysteretic behavior (25|. There- 
fore the current measured in that experiments is proba- 
bly not I c but I r . And indeed, the estimations made in 
Ref. [f| showed that the measured critical current was 
much smaller (more than 10 times) than the depairing 
current. Therefore we believe that our results can be 
directly applied for the explanation of the enhancement 
of the critical current found in Refs. 3-|f|- Taking into 
account that in zinc L, ln ~ 25/_tm [14 1 (when £o — 250 
nm Q) it becomes clear why the effect was weak in rel- 
atively long microbridge with L = 10/im ~ 0ALi n and 
in short one with L — lfj,m ~ 4£o (fj- In Sn the effect 
was absent for nanowires with L = 6 — 35um 0] because 
in tin £o ~ 55 nm and Li n ~ 470 nm [14j . The appli- 
cability of our result to the experiment of Rogachev et. 
al [2J is questionable because those authors claimed that 
their critical current is about the depairing current and 
furthermore experimental IV curves were strongly hys- 
teretic (see inset in Fig. 2(c) in Ref. @). Unfortunately 
no length dependence of the effect was studied in that 
work. 

We believe that our result gives a clue in the under- 
standing of the origin of the negative magnetoresistance 
at low currents / <C I r . It is believed that the finite re- 
sistance of the superconducting nanowires/microbridges 
at low currents originates from the finite rate of thermo- 
activated and/or quantum phase slips (see for example 
review pfjj). Each phase slip event is connected with 
one oscillation of the magnitude of the order parameter 
which provides the 'heating' of the quasiparticles. There- 



From our theoretical calculations we found that the 
enhancement of I r is absent in the temperature interval 
0.92 < T/T c < 0.99 for f in < 20 for all considered lengths 
of the microbridges L = 7 — 51£o- This is not surprising 
because in this temperature interval the corresponding 
Li n is about the coherence length and the period of os- 
cillations of the order parameter is about Ti n . For these 
conditions the effective 'heating' is rather weak (see [221 ]) 
and hence diffusion of 'hot' quasiparticles does not play 
any role. For fi n = 60 the effect is practically absent at 
T/T c > 0.96 (in this case L in ~ 7.5£ ~ £(T)) and it 
becomes noticeable at T/T c — 0.92 for microbridges with 
length 17£ < L < 31£o- 



One more interesting effect which comes from our cal- 
culations is the nonmonotonous dependence of I r on the 
length of the microbridge at zero magnetic field (when 
T' c = T c ). One can see from comparing the values of 
I r in Figs. 3-7 that I r is minimal for L = 21£o- The 
reason for such a behavior is the following - in a very 
short microbridge the 'heating' of quasiparticles is weak 
and I r ~ I c and I c decreases with decreasing L. In a 
very long microbridge L 3> Li n the 'hot' quasiparticles 
relax on the length scale ~ Li n near the phase slip core 
while in microbridge with length L < Li n such a relax- 
ation is less effective. This results in the following main 
tendency: the retrapping current I r first decreases with 
increasing length of the microbridge, reaches the minimal 
value at L > 4£ (when the influence of the leads becomes 
sufficiently weak) and than increases and saturates at 
L > L m . 
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V. CONCLUSION 

The 'heating' of quasiparticles, which occurs in the 
phase slip core due to oscillations of the order param- 
eter can be reduced by diffusion of the quasiparticles to 
the outside regions. It results in an enhancement of the 
retrapping current when one slightly suppresses the or- 
der parameter in the superconducting leads. The en- 
hancement of I r strongly depends on the length of the 
microbridgc - the effect is absent in short mictobridges 
with length L < 2£(T) and it is weak in relatively long 
samples with length L > L; n . Our results predict also 
a nonmonotonous dependence of the retrapping current 
I r on the length of the microbridge - it is minimal when 
4£(T) L < Li n . We should note that our results cannot 



be obtained in the framework of ordinary [25| or extended 
[l8j time-dependent Ginzburg-Landau equations and one 
need to solve Eqs. (1-3) where nonlocal effects connected 
with a time delay of the response and the diffusion of 
noncquilibrium quasiparticles are taken into account. 
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